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Goals  of  Project 

From  our  proposal,  the  four  focus  areas  are: 

•  The  state  of  upper  atmospheric  drivers  during  one  solar  rotation  will  be  different  than  the  state  of  these 
variables  during  the  next.  In  order  to  produce  usable  forecasts,  an  analysis  of  the  inherent  uncertainties  involved 
in  utilizing  the  previous  solar  rotation  drivers  needs  to  be  conducted.  Probability  distribution  functions  (PDFs) 
need  to  be  created  to  better  specify  the  possible  instantaneous  and  temporal  history  of  the  drivers  and  their 
uncertainties. 

•  Because  there  are  large  uncertainties  in  the  predicted  drivers  of  the  system,  ensemble  forecasts  utilizing 
the  PDFs  can  help  in  specifying  the  possible  states  that  may  occur.  The  optimal  technique  for  constructing 
ensembles  of  drivers  from  the  PDFs  to  be  input  into  the  models  using  data  from  the  previous  solar  rotation  will 
need  to  be  determined.  The  number  of  ensemble  members  required  to  adequately  specify  the  forecast  will  need 
to  be  assessed. 

•  The  variability  that  is  lacking  within  empirical  models  of  the  high  latitude  potential  and  aurora  as  well  as 
solar  irradiance  models  will  have  to  be  considered.  Accounting  for  this  variability  will  make  the  model  more 
dynamic  and  realistic. 

•  The  ability  of  an  ensemble-driven  ionosphere-thermosphere  model  to  “predict”  space  weather  needs  to  be 
proven.  This  will  be  done  using  historical  measurements  from  100s  of  ionospheric  and  thermospheric  measure¬ 
ment  locations. 


Accomplishments 

We  have  written  14  papers  with  support  from  this  research  grant,  and  have  completed  the  research  for  at  least 
two  more.  The  14  papers  are  described  with  their  abstracts  below. 

One  of  the  main  studies  that  we  have  conducted  was  the  creation  of  a  new  solar  wind  velocity  precition 
model  based  on  probability  distribution  functions  (PDFs).  The  first  paper  that  was  published  on  this  model 
is  listed  below  and  the  abstract  is  given.  The  second  paper  is  in  preparation  and  will  be  submitted  within  the 
month.  The  abstract  for  that  paper  is  given  here: 

Abrupt  transitions  from  slow  to  fast  solar  wind  represent  a  great  concern  for  the  space  weather  forecasting 
community.  They  cause  important  geomagnetic  storms  that  can  eventually  affect  systems  in  orbit  and  on  the 
ground.  Therefore,  the  Probability  Distribution  Function  (PDF)  model  has  been  improved  to  predict  high 
enhancements  in  the  solar  wind  speed.  New  probability  distribution  functions  allow  for  the  predictions  of  the 
amplitudes  and  the  duration  of  peaks  while  providing  an  interval  of  uncertainty  on  the  prediction.  It  was  found 
that  60%  (91%)  of  the  positive  (negative)  predictions  are  correct.  20%  to  33%  of  the  peaks  in  the  speed  are 
found  by  the  model.  This  represents  a  considerable  improvement  of  the  first  version  of  the  PDF  model.  A  direct 
comparison  to  the  WSA  model  shows  that  the  PDF  model  leads  to  less  false  positive  predictions,  and  miss  less 
events  especially  when  the  peak  reaches  very  high  speeds. 
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Basically,  the  PDF  model  works  in  the  following  way: 

1 .  Given  the  current  solar  wind  speed  and  whether  the  speed  is  increasing  or  decreasing,  the  model  takes 
the  most  probable  speeds  over  the  following  five  days. 

2.  Using  data  from  the  last  day,  the  optimum  delay  between  one  solar  rotation  ago  and  now  can  be  estimated. 
This  only  works  for  a  short  amount  of  time,  though,  since  the  optimal  delay  changes  rapidly. 

3.  The  solar  wind  speed  from  one  solar  rotation  ago  can  be  used  to  predict  the  current  solar  wind  speed 
and  the  solar  wind  speed  for  about  24  hours  into  the  future.  Beyond  that,  the  time  delay  is  most  likely 
incorrect,  and  therefore  the  predictions  based  on  the  one  solar  rotation  ago  values  are  probably  incorrect. 

4.  The  predictions  based  on  one  solar  rotation  ago  and  the  current  speed  and  slope  of  the  speed  are  combined 
linearly  to  allow  for  the  best  possible  solution.  The  coefficients  for  the  two  predictions  are  determined  by 
reducing  the  error  as  compared  to  data. 

5.  The  method  above  is  described  as  the  PDF  model,  and  is  better  than  persistence  after  a  few  hours,  and  is 
statistical!  better  than  the  Wang-Sheely-Arge  (WSA)  model  of  the  solar  wind  velocity. 

6.  The  baseline  PDF  model  does  not  do  a  very  good  job  at  predicting  solar  wind  peaks,  which  is  problematic. 
Therefore,  we  created  an  extension  of  the  PDF  model  that  specifically  looks  for  peaks.  If  fhe  slope  of  fhe 
solar  wind  velocity  exceeds  4.77  km/s/hour,  fhen  fhe  peak-seeking  algorifhm  is  implemenfed. 

7.  There  are  fwo  key  fealures  fhe  peak:  (1)  fhe  lime  of  fhe  peak  and  (2)  fhe  slrenglh  of  fhe  peak.  Probabilily 
dislribulion  functions  are  used  fo  find  each.  The  time  of  fhe  peak  is  eslimafed  by  looking  af  fhe  median 
fime  from  fhe  lowesf  solar  wind  speed  unfil  fhe  peak.  This  time  is  used  in  fhe  model.  Then,  fhe  slrenglh 
of  fhe  peak  is  eslimafed  from  dislribulion  functions  of  Ihe  difference  belween  fhe  peak  speed  and  Ihe 
minimum  in  speed. 

8.  The  currenl  speed  (which  is  Ihe  poinl  al  which  Ihe  slope  exceeded  Ihe  4.77  km/s/hour)  and  Ihe  peak  speed 
are  joined  logelher  wilh  a  polynomial  lhal  is  somewhal  complex  lo  explain  in  a  simple  way. 

9.  If  Ihere  is  a  peak  in  Ihe  magnetic  field  or  Ihe  solar  wind  density  around  Ihe  same  time  as  Ihe  increased 
slope  of  Ihe  solar  wind  velocity,  Ihe  model  does  heller  al  predicting  Ihe  peak  in  velocity,  since  Ihe  increase 
in  Ihe  density  and  magnetic  field  typically  indicates  Ihe  presence  of  a  shock. 

10.  The  ability  of  Ihe  model  lo  predicl  peaks  in  speeds  was  tested  and  compared  lo  models  such  as  Ihe  WSA. 
Il  is  shown  lhal  Ihe  PDF  model  wilh  Ihe  peak  finding  algorilhm  does  belter  lhan  Ihe  WSA  model  Ihe 
majority  of  Ihe  time. 

There  was  significanlly  more  accomplished  on  Ihis  granl,  so  we  lisl  Ihe  papers  published  and  Iheir  abslracls 
below.  This  should  provide  an  overview  of  Ihe  research  lhal  was  conducted. 


Publications  and  Presentations 

Here  is  a  lisl  of  publications  lhal  have  resulted  in  pari  lo  Ibis  granl  (sludenls  are  underlined,  while  posl  docs  are 
in  italics): 
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1.  Liu,  X.,  Ridley,  A.  (2015),  A  simulation  study  of  the  thermosphere  mass  density  response  to  substorms 
using  GITM,  /  Geophys.  Res.  Space  Physics,  120:  7987-8001,  doi:  10.1002/2014JA020962. 

Abstraet:  The  temporal  and  spatial  variations  of  the  thermospherie  mass  density  during  a  series  of  ideal¬ 
ized  substorms  were  investigated  using  the  Global  Ionosphere  Thermosphere  Model  (GITM).  The  max¬ 
imum  mass  density  perturbation  of  an  idealized  substorm  with  a  peak  variation  of  hemispherie  power 
(HP)  of  50  GW  and  interplanetary  magnetie  field  (IMF)  Bz  of  -2  nT  was  14%  about  50  min  after  the 
substorm  onset  in  the  nightside  seetor  of  the  auroral  zone.  The  mass  density  response  to  different  types 
of  energy  input  has  a  strong  loeal  time  dependenee,  with  the  mass  density  perturbation  due  to  only  an 
IMF  Bz  variation  peaking  in  the  dusk  seetor  and  the  density  perturbation  due  to  only  HP  variation  peaks 
in  the  nightside  seetor.  Simulations  with  IMF  Bz  ehanges  only  and  HP  ehanges  only  showed  that  the 
system  behaves  slightly  nonlinearly  when  both  IMF  and  HP  variations  are  ineluded  (a  maximum  of  6% 
of  the  nonlinearity)  and  that  the  nonlinearity  grows  with  energy  input.  The  neutral  gas  heating  rate  due 
to  Joule  heating  was  of  same  magnitude  as  the  heating  rate  due  to  preeipitation,  but  the  majority  of  the 
temperature  enhaneement  due  to  the  heating  due  to  preeipitation  oeeurs  at  lower  altitude  as  eompared  to 
the  auroral  heating.  About  110  min  after  onset,  a  negative  mass  density  perturbation  ( -5%)  oeeurred  in 
the  night  seetor,  whieh  was  eonsistent  with  the  mass  density  measurement  of  the  CHAMP  satellite. 

2.  Wang,  H.,  Ridley,  A.,  Zhu,  J.  (2015),  Theoretieal  study  of  zonal  differenees  of  eleetron  density 
at  midlatitudes  with  GITM  simulation,  J.  Geophys.  Res.  Space  Physics,  120,  2951-2966,  doi: 
10.1002/2014JA020790. 

Abstraet:  This  study  investigated  various  physieal  proeesses  responsible  for  the  longitudinal  modula¬ 
tion  of  eleetron  density  (Ne)  at  midlatitudes  by  employing  the  global  ionosphere-thermosphere  model 
(GITM).  The  good  agreements  between  GITM  outputs  and  CHAMP  observations  indieate  that  the  model 
is  a  suitable  tool  to  perform  the  theoretieal  study.  Nine  runs  were  carried  out  to  determine  the  effects 
from  geomagnetic  field  geometry,  zonal  wind,  meridional  wind,  high-latitude  activity,  migrating  tides 
from  the  lower  atmosphere,  and  solar  illumination  in  quantitative  ways.  Distinct  features  were  discussed 
as  follows.  It  was  crucial  that  the  geomagnetic  and  geographical  axes  were  offset  for  the  development 
of  the  longitudinal  difference  of  Ne.  The  zonal  wind  contributes  to  about  80%  of  the  fraction  of  the 
observed  longitudinal  dependence  of  Ne.  The  meridional  wind  effect  is  out  of  phase  with  the  zonal  wind 
over  North  America  and  Southern  Ocean  regions,  which  trims  the  fraction  of  the  longitudinal  difference 
to  65%.  Over  the  South  Pacific  Ocean,  the  nighttime  Ne  maintains  at  a  higher  level  because  of  in-phase 
effects  from  both  zonal  and  meridional  winds.  The  solar  illumination  was  important  in  the  formation 
of  the  background  longitudinal  pattern  of  the  electron  density.  The  migrating  tide  from  the  lower  atmo¬ 
sphere  could  enhance  the  longitudinal  difference  of  Ne  by  15%  over  North  America.  Enhanced  activities 
at  high  latitudes  could  alter  the  longitudinal  pattern  of  Ne  by  transporting  thermospheric  composition 
disturbances  to  midlatitudes. 

3.  Sheng,  C.,  Deng,  Y,  Wu,  Q.,  Ridley,  A.  and  Haggstrdm,  I.  (2015),  Thermospheric  winds  around  the  cusp 
region.  J.  Geophys.  Res.  Space  Physics,  120:  12481255.  doi:  10.1002/2014JA020028. 

Abstract:  An  equatorward  wind  has  been  observed  first  by  the  balloon-borne  Fabry-Perot  interferometer 
called  High- Altitude  Interferometer  Wind  Observation  on  the  equatorward  side  of  the  cusp  near  the  local 
noon,  which  is  opposite  to  the  typical  direction  of  neutral  wind  driven  by  the  day-night  pressure  gradient. 
However,  this  dayside  equatorward  wind  was  not  reproduced  by  the  standard  Thermosphere  Ionosphere 
Electrodynamics  General  Circulation  Model  under  the  resolution  of  5°  longitude  by  5°  latitude  (5°  x 
5°).  In  this  study,  the  Global  Ionosphere  Thermosphere  Model  has  been  run  in  different  cases  and  un¬ 
der  different  resolutions  to  investigate  the  neutral  dynamics  around  the  cusp  region.  Eirst,  we  compare 
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the  simulations  with  and  without  additional  cusp  energy  inputs  to  identify  the  influence  of  cusp  heating. 
Both  runs  have  a  resolution  of  5°  x  1°  (longitude  latitude)  in  order  to  better  resolve  the  cusp  region.  Af¬ 
ter  adding  in  the  cusp  energy,  the  meridional  wind  in  simulation  turns  to  be  equatorward  on  the  dayside, 
which  is  consistent  with  the  observation.  It  indicates  that  strong  heating  in  the  cusp  region  causes  changes 
in  the  pressure  gradient  around  the  cusp  and  subsequent  variations  in  the  neutral  winds.  The  simulations 
with  the  same  cusp  heating  specifications  are  repeated,  but  with  different  horizontal  resolutions  to  exam¬ 
ine  the  influence  of  resolution  on  the  simulation  results.  The  comparisons  show  that  the  resolution  of  5° 
X  1°  can  resolve  the  cusp  region  much  more  stably  and  consistently  than  the  5°  x  5°  resolution. 

4.  Mannucci,  A.  J.,  O.  P.  Verkhoglyadova,  B.  T.  Tsurutani,  X.  Meng,  X.  Pi,  C.  Wang,  G.  Rosen,  E.  Lynch, 

S.  Sharma,  A.  Ridley,  W.  Manchester,  B.  Van  Der  Holst,  E.  Echer,  and  R.  Hajra  (2015),  Medium-Range 
Thermosphere-Ionosphere  Storm  Eorecasts.  Space  Weather,  13,  125129.  doi:  10.1002/2014SW001125. 

5.  Zou,  S.,  M.  B.  Moldwin,  A.  J.  Ridley,  M.  J.  Nicolls,  A.  J.  Coster,  E.  G.  Thomas,  and  J.  M.  Ruohoniemi 
(2014),  On  the  generation/decay  of  the  storm-enhanced  density  plumes:  Role  of  the  convection  flow  and 
field-aligned  ion  flow,  7.  Geophys.  Res.  Space  Physics,  119,  85438559,  doi:  10. 1002/2014JA020408. 

Abstract:  Storm-enhanced  density  (SED)  plumes  are  prominent  ionospheric  electron  density  increases  at 
the  dayside  middle  and  high  latitudes.  The  generation  and  decay  mechanisms  of  the  plumes  are  still  not 
clear.  We  present  observations  of  SED  plumes  during  six  storms  between  2010  and  2013  and  compre¬ 
hensively  analyze  the  associated  ionospheric  parameters  within  the  plumes,  including  vertical  ion  flow, 
field-aligned  ion  flow  and  flux,  plasma  temperature,  and  field-aligned  currents,  obtained  from  multiple 
instruments,  including  GPS  total  electron  content  (TEC),  Poker  Plat  Incoherent  Scatter  Radar  (PPISR), 
Super  Dual  Auroral  Radar  Network,  and  Active  Magnetosphere  and  Planetary  Electrodynamics  Response 
Experiment.  The  TEC  increase  within  the  SED  plumes  at  the  PPISR  site  can  be  1.4-5. 5  times  their  quiet 
time  value.  The  plumes  are  usually  associated  with  northwestward  ExB  flows  ranging  from  a  couple  of 
hundred  m  s-1  to  >  1  km  s-1.  Upward  vertical  flows  due  to  the  projection  of  these  ExB  drifts  are  mainly 
responsible  for  lifting  the  plasma  in  sunlit  regions  to  higher  altitude  and  thus  leading  to  plume  density 
enhancement.  The  upward  vertical  flows  near  the  poleward  part  of  the  plumes  are  more  persistent,  while 
those  near  the  equatorward  part  are  more  patchy.  In  addition,  the  plumes  can  be  collocated  with  either 
upward  or  downward  field-aligned  currents  (PACs)  but  are  usually  observed  equatorward  of  the  peak 
of  the  Region  1  upward  PAG,  suggesting  that  the  northwestward  flows  collocated  with  plumes  can  be 
either  subauroral  or  auroral  flows.  Purthermore,  during  the  decay  phase  of  the  plume,  large  downward 
ion  flows,  as  large  as  approximately  200  m  s-1,  and  downward  fluxes,  as  large  as  10^^  m-2  s-1,  are  often 
observed  within  the  plumes.  In  our  study  of  six  storms,  enhanced  ambipolar  diffusion  due  to  an  elevated 
pressure  gradient  is  able  to  explain  two  of  the  four  large  downward  flow/flux  cases,  but  this  mechanism 
is  not  sufficient  for  the  other  two  cases  where  the  flows  are  of  larger  magnitude.  Por  the  latter  two  cases, 
enhanced  poleward  thermospheric  wind  is  suggested  to  be  another  mechanism  for  pushing  the  plasma 
downward  along  the  field  line.  These  downward  flows  should  be  an  important  mechanism  for  the  decay 
of  the  SED  plumes. 

6.  Bussy-Virat,  C.  D.,  and  A.  J.  Ridley,  Predictions  of  the  solar  wind  speed  by  the  probability  distribution 
function  model.  Space  Weather,  12,  337353,  doi:10.1002/2014SW001051,  2014. 

Abstract:  The  near-Earth  space  environment  is  strongly  driven  by  the  solar  wind  and  interplanetary  mag¬ 
netic  field.  This  study  presents  a  model  for  predicting  the  solar  wind  speed  up  to  5  days  in  advance. 
Probability  distribution  functions  (PDPs)  were  created  that  relate  the  current  solar  wind  speed  and  slope 
to  the  future  solar  wind  speed,  as  well  as  the  solar  wind  speed  to  the  solar  wind  speed  one  solar  rota¬ 
tion  in  the  future.  It  was  found  that  a  major  limitation  of  this  type  of  technique  is  that  the  solar  wind 
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periodicity  is  close  to  27  days  but  can  be  from  about  22  to  32  days.  Further,  the  optimum  lag  between 
two  solar  rotations  can  change  from  day  to  day,  making  a  prediction  of  the  future  solar  wind  speed  based 
solely  on  the  solar  wind  speed  approximately  27  days  ago  quite  difficult.  It  was  found  that  using  a  linear 
combination  of  the  solar  wind  speed  one  solar  rotation  ago  and  a  prediction  of  the  solar  wind  speed  based 
on  the  current  speed  and  slope  is  optimal.  The  linear  weights  change  as  a  function  of  the  prediction 
horizon,  with  shorter  prediction  times  putting  more  weight  on  the  prediction  based  on  the  current  solar 
wind  speed  and  the  longer  prediction  times  based  on  an  even  spread  between  the  two.  For  all  prediction 
horizons  from  8  h  up  to  120  h,  the  PDF  Model  is  shown  to  be  better  than  using  the  current  solar  wind 
speed  (i.e.,  persistence),  and  better  than  the  Wang-Sheeley-Arge  Model  for  prediction  horizons  of  24  h. 

7.  Deng,  Y.,  and  A.  J.  Ridley,  Simulation  of  non-hydrostatic  gravity  wave  propagation  in  the  upper  atmo¬ 
sphere,  Annales  Geophysicae,  32,  443447,  doi:10.5194/angeo-32-443-2014,  2014. 

Abstract:  The  high-frequency  and  small  horizontal  scale  gravity  waves  may  be  reflected  and  ducted  in 
non-hydrostatic  simulations,  but  usually  propagate  vertically  in  hydrostatic  models.  To  examine  gravity 
wave  propagation,  a  preliminary  study  has  been  conducted  with  a  global  ionosphere-thermosphere  model 
(GITM),  which  is  a  non-hydrostatic  general  circulation  model  for  the  upper  atmosphere.  GITM  has 
been  run  regionally  with  a  horizontal  resolution  of  0.2°  long  0.2°  lat  to  resolve  the  gravity  wave  with 
wavelength  of  250  km.  A  cosine  wave  oscillation  with  amplitude  of  30  m  s-1  has  been  applied  to  the 
zonal  wind  at  the  low  boundary,  and  both  high-frequency  and  low-frequency  waves  have  been  tested.  In 
the  high-frequency  case,  the  gravity  wave  stays  below  200  km,  which  indicates  that  the  wave  is  reflected 
or  ducted  in  propagation.  The  results  are  consistent  with  the  theoretical  analysis  from  the  dispersion 
relationship  when  the  wavelength  is  larger  than  the  cutoff  wavelength  for  the  non-hydrostatic  situation. 
However,  the  low-frequency  wave  propagates  to  the  high  altitudes  during  the  whole  simulation  period, 
and  the  amplitude  increases  with  height.  This  study  shows  that  the  non-hydrostatic  model  successfully 
reproduces  the  high-frequency  gravity  wave  dissipation. 

8.  Ridley,  A.  J.,  A.  M.  Dodger,  and  M.  W.  Liemohn,  Exploring  the  efficacy  of  different  electric  field  models 
in  driving  a  model  of  the  plasmasphere,  J.  Geophys.  Res.,  119,  46214638,  doi:  10.1002/2014JA019836, 
2014. 

Abstract:  The  dynamics  of  the  plasmasphere  are  strongly  controlled  by  the  inner  magnetospheric  electric 
field.  In  order  to  capture  realistically  the  erosion  of  the  nightside  plasmapause  and  the  formation  of  the 
drainage  plume  in  a  model  of  the  plasmasphere,  the  electric  field  must  be  accurate.  This  study  investi¬ 
gates  how  well  five  different  electric  field  models  drive  the  Dynamic  Global  Core  Plasma  Model  during 
eight  storm  periods.  The  five  electric  field  models  are  the  Volland-Stern  analytic  formula  with  Maynard- 
Chen  Kpdependence,  two  versions  of  the  Weimer  statistical  models  (96  and  05),  and  two  versions  of 
the  Assimilative  Mapping  of  Ionospheric  Electrodynamics  (AMIE)  technique  using  magnetometer  and 
DMSP  satellite  data.  Manually  extracted  plasmapause  locations  from  images  taken  by  the  EUV  instru¬ 
ment  on  the  Imager  for  Magnetopause-to-Aurora  Global  Exploration  (IMAGE)  satellite,  as  described  by 
Goldstein  et  al.  (2005),  were  compared  to  the  simulation  results  throughout  the  main  phase  of  the  eight 
events.  Three  methods  of  calculating  the  plasmapause  were  employed  to  determine  the  best  fit  to  EUV 
data,  using  the  maximum  gradient,  a  constant  density  contour  (fit  method),  and  the  location  in  which 
the  modeled  density  fell  significantly  below  the  specified  saturation  density  for  the  given  radial  position 
(saturation  method).  It  was  found  that  the  simulations  driven  by  the  Weimer  (1996)  model  produced  the 
best  fit  overall  and  that  the  fit  and  saturation  methods  worked  best  for  matching  the  model  results  to  the 
observations. 

9.  Zou,  S.,  A.  J.  Ridley,  M.  B.  Moldwin,  M.  J.  Nicolls,  A.  J.  Coster,  E.  G.  Thomas,  and  J.  M.  Ruohoniemi, 
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Multi-instrument  observations  of  SED  during  2425  Oetober  2011  storm:  Implieations  for  SED  formation 
proeesses,  J.  Geophys.  Res.,  118,  77987809,  doi:10.1002/2013JA018860,  2013. 

Abstraet:present  multiple  instrument  observations  of  a  storm-enhaneed  density  (SED)  during  the  24-25 
October  2011  intense  geomagnetic  storm.  Eormation  and  the  subsequent  evolution  of  the  SED  and  the 
midlatitude  trough  are  revealed  by  global  GPS  vertical  total  electron  content  maps.  In  addition,  we 
present  high  time  resolution  Poker  Plat  Incoherent  Scatter  Radar  (PPISR)  observations  of  ionospheric 
profiles  within  the  SED.  We  divided  the  SED  observed  by  PPISR  into  two  parts.  Both  parts  are  charac¬ 
terized  by  elevated  ionospheric  peak  height  (hmP2)  and  total  electron  content,  compared  to  quiet  time 
values.  However,  the  two  parts  of  the  SED  have  different  characteristics  in  the  electron  temperature  (Te), 
the  P  region  peak  density  (NmP2),  and  convection  flows.  The  first  part  of  the  SED  is  associated  with 
enhanced  Te  in  the  lower  P  region  and  reduced  Te  in  the  upper  P  region  and  is  collocated  with  northward 
convection  flows.  The  NmP2  was  lower  than  quiet  time  values.  The  second  part  of  the  SED  is  associ¬ 
ated  with  significantly  increased  NmP2,  elevated  Te  at  all  altitudes  and  is  located  near  the  equatorward 
boundary  of  large  northwestward  flows.  Based  on  these  observations,  we  suggest  that  the  mechanisms 
responsible  for  the  formation  of  the  two  parts  of  the  SED  may  be  different.  The  first  part  is  due  to 
equatorward  expansion  of  the  convection  pattern  and  the  projection  of  northward  convection  flows  in  the 
vertical  direction,  which  lifts  the  ionospheric  plasma  to  higher  altitudes  and  thus  reduces  the  loss  rate 
of  plasma  recombination.  The  second  part  is  more  complicated.  Besides  equatorward  expansion  of  the 
convection  pattern  and  large  upward  flows,  evidences  of  other  mechanisms,  including  horizontal  advec- 
tion  due  to  fast  flows,  energetic  particle  precipitation,  and  enhanced  thermospheric  wind  in  the  topside 
ionosphere,  are  also  present.  Estimates  show  that  contribution  from  precipitating  energetic  protons  is  at 
most  approximately  10%  of  the  total  E  region  density.  The  thermospheric  wind  also  plays  a  minor  role 
in  this  case. 

10.  Pulkkinen,  A.,  E.  Rastatter,  M.  Kuznetsova,  H.  Singer,  C.  Balch,  D.  Weimer,  G.  Toth,  A.  Ridley,  T. 
Gombosi,  M.  Wiltberger,  J.  Raeder,  R.  Weigel,  Community-wide  validation  of  geospace  model  ground 
magnetic  field  perturbation  predictions  to  support  model  transition  to  operations.  Space  Weather,  11, 
369-385,  doi:  10.1002/swe.20056,  2013. 

Abstract:  In  this  paper  we  continue  the  community-wide  rigorous  modern  space  weather  model  validation 
efforts  carried  out  within  GEM,  CEDAR  and  SHINE  programs.  In  this  particular  effort,  in  coordination 
among  the  Community  Coordinated  Modeling  Center  (CCMC),  NCAA  Space  Weather  Prediction  Center 
(SWPC),  modelers,  and  science  community,  we  focus  on  studying  the  models’  capability  to  reproduce 
observed  ground  magnetic  field  fluctuations,  which  are  closely  related  to  geomagnetically  induced  cur¬ 
rent  phenomenon.  One  of  the  primary  motivations  of  the  work  is  to  support  NOAA  SWPC  in  their 
selection  of  the  next  numerical  model  that  will  be  transitioned  into  operations.  Six  geomagnetic  events 
and  12  geomagnetic  observatories  were  selected  for  validation.  While  modeled  and  observed  magnetic 
field  time  series  are  available  for  all  12  stations,  the  primary  metrics  analysis  is  based  on  six  stations 
that  were  selected  to  represent  the  high-latitude  and  mid-latitude  locations.  Events-based  analysis  and 
the  corresponding  contingency  tables  were  built  for  each  event  and  each  station.  The  elements  in  the 
contingency  table  were  then  used  to  calculate  Probability  of  Detection  (POD),  Probability  of  Ealse  De¬ 
tection  (POED)  and  Heidke  Skill  Score  (HSS)  for  rigorous  quantification  of  the  models’  performance.  In 
this  paper  the  summary  results  of  the  metrics  analyses  are  reported  in  terms  of  POD,  POED  and  HSS. 
More  detailed  analyses  can  be  carried  out  using  the  event  by  event  contingency  tables  provided  as  an 
online  appendix.  An  online  interface  built  at  CCMC  and  described  in  the  supporting  information  is  also 
available  for  more  detailed  time  series  analyses. 

11.  Honkonen,  I.,  E.  Rastatter,  A.  Grocott,  A.  Pulkkinen,  M.  Palmroth,  J.  Raeder,  A.J.  Ridley,  M.  Wilt- 
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berger,  On  the  performance  of  global  magnetohydrodynamic  models  in  the  Earth’s  magnetosphere,  Space 
Weather,  11,  313-326,  doi:  10.1002/swe.20055,  2013. 

Abstract:  We  study  the  performance  of  four  magnetohydrodynamic  models  (BATS-R-US,  GUMICS, 
LEM,  OpenGGCM)  in  the  Earth’s  magnetosphere.  Using  the  Community  Coordinated  Modeling  Cen¬ 
ter’s  Run-on-Request  system,  we  compare  model  predictions  with  magnetic  field  measurements  of  the 
Cluster,  Geotail  and  Wind  spacecraft  during  a  multiple  substorm  event.  We  also  compare  model  cross 
polar  cap  potential  results  to  those  obtained  from  the  Super  Dual  Auroral  Radar  Network  (SuperDARN) 
and  the  model  magnetopause  standoff  distances  to  an  empirical  magnetopause  model.  The  correlation 
coefficient  (CC)  and  prediction  efficiency  (PE)  mefrics  are  used  fo  objectively  evaluafe  model  perfor¬ 
mance  quanfifafively.  Eor  all  four  models,  fhe  besf  performance  oufside  geosynchronous  orbif  is  found 
on  fhe  dayside.  Generally,  fhe  performance  of  models  decreases  sfeadily  downsfream  from  fhe  Earfh.  On 
fhe  dayside  mosf  CCs  are  above  0.5  wifh  CCs  for  Bx  and  Bz  close  fo  0.9  for  fhree  ouf  of  four  models.  In 
fhe  magnefofail  af  a  disfance  of  abouf  -130  Earfh  radii  from  Earfh,  fhe  predicfion  efficiency  of  all  models 
is  below  fhaf  of  using  an  average  value  for  fhe  prediction  wifh  fhe  excepfion  of  Bz.  Bx  is  mosf  oflen  besf 
predicted  and  correlated  bofh  on  fhe  dayside  and  fhe  nighfside  close  fo  fhe  Earfh  whereas  in  fhe  far  fail 
fhe  CC  and  PE  for  Bz  are  subsfanfially  higher  fhan  ofher  componenfs  in  all  models.  We  also  find  fhaf 
increasing  fhe  resolufion  or  coupling  an  additional  physics  module  does  nol  aufomafically  increase  fhe 
model  performance  in  fhe  magnefosphere. 

12.  Wei,  Y.,  W.  Wan,  B.  Zhao,  M.  Hong,  A.  Ridley,  Z.  Ren,  M.  Eraenz,  E.  Dubinin,  and  M.  He,  Solar  wind 
densify  confrolling  penefrafion  elecfric  field  af  fhe  equaforial  ionosphere  during  a  safurafion  of  cross  polar 
cap  potential,  J.  Geophys.  Res.,  117,  A09308,  doi:10.1029/2012JA017597,  2012. 

Abslracf:  The  mosf  imporfanf  source  of  elecfrodynamic  disfurbances  in  fhe  equatorial  ionosphere  during 
fhe  main  phase  of  a  sform  is  fhe  prompf  penefrafion  elecfric  field  (PPEE)  originating  from  fhe  high-lafifude 
region.  If  has  been  known  fhaf  such  an  elecfric  field  is  correlafed  wifh  fhe  magnefospheric  convecfion  or 
inferplanefary  elecfric  field.  Here  we  show  a  unique  case,  in  which  fhe  elecfric  field  disfurbance  in  fhe 
equaforial  ionosphere  cannof  be  interpreted  by  fhis  concepf.  During  fhe  supersform  on  Nov.  20-21,  2003, 
fhe  cross  polar  cap  pofenfial  (CPCP)  safurafed  af  leasf  for  8.2  h.  The  CPCP  reconsfrucfed  by  Assimilative 
Mapping  of  Ionospheric  Electrodynamics  (AMIE)  procedure  suggested  that  the  PPEE  at  the  equatorial 
ionosphere  still  correlated  with  the  saturated  CPCP,  but  the  CPCP  was  controlled  by  the  solar  wind 
density  instead  of  the  interplanetary  electric  field.  However,  fhe  predicfed  CPCPs  by  Hill-Siscoe-Ober 
(HSO)  model  and  Boyle-Ridley  (BR)  model  were  nol  fully  consisfenf  wifh  fhe  AMIE  resulf  and  PPEE. 
The  PPEE  also  decoupled  from  fhe  convecfion  elecfric  field  in  fhe  magnefolail.  Due  fo  fhe  decoupling, 
fhe  elecfric  field  in  fhe  ring  currenf  was  nol  able  fo  comply  wifh  fhe  varialions  of  PPEE,  and  fhis  resulled 
in  a  long-duralion  elecfric  field  penefrafion  wilhoul  shielding. 

13.  J.S.  Shim,  Kuznelsova,  M.;  Raslfler,  E.;  Bililza,  D.;  Bulala,  M.;  Codrescu,  M.;  Emery,  B.  A.;  Eosler, 
B.;  Euller-Rowell,  T.  J.;  Huba,  J.;  Mannucci,  A.  J.;  Pi,  X.;  Ridley,  A.;  Scherliess,  E.;  Schunk,  R.  W.; 
Sojka,  J.  J.;  Stephens,  P;  Thompson,  D.  C.;  Weimer,  D.;  Zhu,  E.;  Sulfon,  E.,  CEDAR  Eleclrodynamics 
Thermosphere  Ionosphere  (ETI)  Challenge  for  syslemalic  assessmenl  of  ionosphere/lhermosphere  mod¬ 
els:  Eleclron  densify,  neulral  densify,  NmE2,  and  hmE2  using  space  based  observalions.  Space  Weather, 
10,  doi:  10.1 029/20 12SW000851,  2012. 

Abslracf:  In  an  efforl  fo  quanfifafively  assess  fhe  currenf  capabilities  of  lonosphere/Thermosphere  (IT) 
models,  an  IT  model  validation  sludy  using  mefrics  was  performed.  This  sludy  is  a  main  pari  of  fhe 
CEDAR  Eleclrodynamics  Thermosphere  Ionosphere  (ETI)  Challenge,  which  was  inilialed  af  fhe  CEDAR 
workshop  in  2009  fo  heifer  comprehend  slrenglhs  and  weaknesses  of  models  in  predicting  fhe  IT  syslem. 
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and  to  trace  improvements  in  ionospheric/thermospheric  specification  and  forecast.  For  the  challenge, 
two  strong  geomagnetic  storms,  four  moderate  storms,  and  three  quiet  time  intervals  were  selected.  For 
the  selected  events,  we  obtained  four  scores  (i.e.,  RMS  error,  prediction  efficiency,  rafio  of  fhe  maximum 
change  in  amplifudes,  and  rafio  of  fhe  maximum  amplifudes)  fo  compare  fhe  performance  of  models 
in  reproducing  fhe  selecfed  physical  parameters  such  as  verfical  drifls,  elecfron  and  neufral  densifies, 
NmF2,  and  hmF2.  In  fhis  paper,  we  presenf  fhe  resulfs  from  comparing  modeled  values  againsf  space- 
based  measuremenfs  including  NmF2  and  hmF2  from  fhe  CHAMP  and  COSMIC  safelliles,  and  elecfron 
and  neufral  densities  af  fhe  CHAMP  safellile  locations.  If  is  found  fhaf  fhe  accuracy  of  models  varies  wifh 
fhe  mefrics  used,  lafifude  and  geomagnetic  acfivify  level. 

14.  Yigit,  E.,  A.  J.  Ridley,  and  M.  B.  Moldwin,  Imporfance  of  capfuring  heliospheric  variabilify  for  sfudies 
of  fhermospheric  vertical  winds,  J.  Geophys.  Res.,  117,  A07306,  doi:10.1029/2012JA017596,  2012. 

Abslracf:  Using  fhe  Global  Ionosphere  Thermosphere  Model  wifh  observed  real-lime  heliospheric  inpuf 
dafa,  fhe  magnifude  and  variabilify  of  fhermospheric  neufral  vertical  winds  are  invesfigafed.  In  order 
fo  determine  fhe  role  of  variabilify  in  fhe  Inlerplanelary  Magnetic  Field  (IMF)  and  solar  wind  densify 
on  fhe  neufral  wind  variabilify,  fhe  heliospheric  inpuf  dafa  are  smoolhed.  The  effecfs  of  smoolhing  fhe 
IMF  and  solar  wind  and  densify  on  fhe  vertical  winds  are  simulated  for  fhe  cases  of  no  smoolhing,  5- 
minule,  and  12-minule  smoolhing.  Various  verfical  wind  acceleration  lerms,  such  as  fhe  nonhydroslalic 
acceleralion,  are  quanlified.  Polar  slereographic  projeclions  of  fhe  variabililies  of  verfical  wind  and  ion 
flows  are  compared  fo  highlighl  exisling  correlations.  Overall,  fhe  smoolher,  lhal  is,  fhe  less  variable  fhe 
IMF  and  solar  wind  paramelers  are,  fhe  weaker  are  fhe  magnifude  and  fhe  variabilify  of  fhe  fhermospheric 
vertical  winds.  Weaker  IMF  variabilify  leads  fo  smaller  variabilify  in  ion  flows,  which  in  lurn  negatively 
impacls  fhe  variabilify  and  fhe  magnifude  of  Joule  heating.  Small-scale  lemporal  varialion  of  fhe  verfical 
wind  acceleralion,  and  Ihus  fhe  variabilify  of  fhe  verfical  wind,  is  dominated  by  fhe  nonhydroslalic  lerm 
lhal  is  conlrolled  primarily  by  fhe  lemporal  variation  of  fhe  Joule  heating,  which  in  lurn  is  relaled  fo  ion 
flow  variations  lhal  are  shaped  by  fhe  IMF  in  fhe  high-lalilude  Ihermosphere.  Wavelel  analysis  of  fhe 
vertical  wind  dafa  shows  lhal  gravity  waves  of  5  and  10-minule  periods  are  more  prominenl  when  fhe 
model  is  run  wifh  high-resolulion  real-time  IMF  and  solar  wind  dafa.  Belter  capfuring  of  fhe  lemporal 
varialion  of  fhe  IMF  and  solar  wind  paramelers  is  crucial  for  modeling  fhe  variabilify  and  magnifude  of 
fhermospheric  verfical  winds. 
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In  addition  to  these  talks,  we  have  given  many  talks  at  meetings  such  as  GEM,  CEDAR,  and  other  workshop- 
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